A self-sustaining nonlinear dynamo process in Keplerian shear flows 
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A three-dimensional nonlinear dynamo process is identified in rotating plane Couette flow in the Keplerian 
regime. It is analogous to the hydrodynamic self-sustaining process in non-rotating shear flows and relies on the 
magneto-rotational instability of a toroidal magnetic field. Steady nonlinear solutions are computed numerically 
for a wide range of magnetic Reynolds numbers but are restricted to low Reynolds numbers. This process may 
be important to explain the sustenance of coherent fields and turbulent motions in Keplerian accretion disks, 
where all its basic ingredients are present. 
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The most natural explanation for the efficient outward 
angular momentum transport inferred from observed ac- 
cretion luminosities in accretion disks is that these objects 
are turbulent In magnetized disks, turbulence is likely 
triggered by the magneto-rotational instability [MRI, see 
When a mean magnetic field with non-zero net flux 
perpendicular to the disk plane is imposed, the instability 
takes the form of a two-dimensional channel flow which 
breaks down into three-dimensional MHD turbulence O] 
as a result of secondary Kelvin-Helmholtz instabilities (j]. 
This simple field configuration, however, may not be rel- 
evant to all accretion disks. Recent observations have for 
instance demonstrated that the magnetic fields in the inner- 
most regions of some disks are probably not created by an 
external object such as the central accreting body but are 
intrinsic to the disk ISj]. Their precise origin, however, re- 
mains largely unknown. Understanding the nonlinear sat- 
uration of the MRI in the absence of an externally imposed 
magnetic field or more generally when there is no net mag- 
netic flux threading through the disk proves to be a com- 
plicated task because one needs to explain in that case how 
the magnetic field, whose presence within the disk is per- 
manently required for the generation of MHD turbulence, 
can be sustained against dissipation. Local numerical disk 
simulations in the absence of a net magnetic flux I^lZJ sug- 
gest that a dynamo process possibly relying on the MRI of 
toroidal magnetic fields fsl] may be at work, but a detailed 
theoretical understanding of such a scenario is still lacking. 

In this Letter, we report the discovery of self-sustaining 
dynamo action in magnetized, spanwise rotating plane 
Couette flow (PCF) in the Keplerian regime (linearly stable 
from the purely hydrodynamic point of view) characteristic 
of accretion disks. The phenomenology of this magneto- 
hydrodynamic (MHD) process is analogous to that of the 
self-sustaining process (SSP) thought to be responsible for 
the transition to turbulence in hydrodynamic shear flows 
dslll^l- It relies on three fundamental physical effects: 

1. linear amplification of zero-net-flux toroidal (az- 
imuthal in the accretion disk terminology) magnetic 
field induced by the distortion of a weak poloidal 



(radial and vertical) seed field by the background 
shear (differential rotation) ; 

2. three-dimensional linear instability (MRI) of the 
toroidal magnetic field ; 

3. regeneration of the poloidal field owing to the non- 
linear feedback of the MRI modes. 

Even though linear processes are essential to the mech- 
anism, either via transient linear growth or linear insta- 
bility, the whole process is fundamentally nonlinear since 
the self-sustaining loop cannot be closed without nonlinear 
feedback. This notably means that there is no kinematic 
regime for the dynamo, which is therefore subcritical: an 
initial finite-amplitude, zero-net-flux magnetic field distur- 
bance is required for the dynamo to operate. This con- 
trasts with the externally imposed field problem in which 
infinitesimal disturbances grow exponentially due to the 
MRI. To the best of our knowledge, this is the first in- 
stance of an explicit nonlinear subcritical dynamo solution. 
We first describe the three steps of the process, construct 
steady nonlinear solutions for Keplerian PCF at various 
magnetic and kinetic Reynolds numbers and finally discuss 
the relevance of our results to accretion disk theory. 

We consider PCF for an incompressible fluid with unit 
density, constant kinematic viscosity v and magnetic diffu- 
sivity rj. The flow is driven by two counter-moving rigid, 
no-slip, perfectly conducting walls located at y = +d and 
is rotating at a constant rate Q. along the spanwise (ver- 
tical) z-axis perpendicular to the linear background shear 
flow Vb (y) - SyCx- It is taken to be spatially periodic in 
both streamwise (toroidal, x) and spanwise (z) directions, 
with periods L j = Ind/a and L, = 27td//3. Using 1/S as 
a time unit and the channel half-width d as a length unit, 
we define a Reynolds number Re - Sd^/v and a mag- 
netic Reynolds number Rm - Sd^/rj. This configuration 
represents an idealized local model of a differentially ro- 
tating Keplerian accretion disk provided that the rotation 
number 7?n = -2Q./S equals -4/3 (anticyclonic Rayleigh- 
stable rotation). Owing to PCF symmetries, we look for 
nonUnear three-dimensional steady solutions of the incom- 
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pressible MHD equations (Navier-Stokes equation with a 
Lorentz force and induction equation) for magnetic and ve- 
locity field perturbations b and v of Keplerian PCF. We first 
consider step 1 of the process. We define a streamfunction 
iff(y, z, t) and a flux function xiy, z, t) to describe the mean 
poloidal fields Vp = V x (_\]/ex) and fip = V x ix^x), where 
overbars stand for x-averaging, and denote the x-dependent 
part of the fields by v' and b' . The toroidal and poloidal 
components of the x-averaged induction equation read 
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Equation ([T]) has a linear induction source term by due to 
the presence of the background shear, which makes it pos- 
sible to generate an (9(1) toroidal magnetic field from a 
weak (9(1 /Rm) poloidal field (this results mathematically 
from the non-normality of the linear operator). This ax- 
isymmetric process for two-dimensional, three-component 
(2D-3C) fields is often referred to as the Q-efFect in dy- 
namo theory and is the MHD analog of the algebraic am- 
plification of streaks (streamwise velocity field) in non- 
rotating hydrodynamic shear flows, known as the lift-up 
eff'ect in the fluid dynamics literature 111 ill . Unless some 
three-dimensional nonlinear mechanism regenerates the 
poloidal field, this linear process can however only be tran- 
sient: for purely two-dimensional configurations, Eq. (|2]l 
is a simple advection-diffusion equation with no source 
term for x, which must therefore decay resistively on an 
(9(Rm) timescale. This is just a restatement of Cowl- 
ing's theorem for an axisymmetric (x-independent here) 
system. To obtain a SSP, we must therefore consider the 
x-dependent, three-dimensional instabilities of the tran- 
siently amplified toroidal component of the 2D-3C field 
(step 2) and their nonlinear feedback in the poloidal equa- 
tion (step 3). Three-dimensional, nonlinear steady solu- 
tions are possible only if the nonlinear interaction term on 
the r.h.s. of Eq. (|2]i has the ability to regenerate the original 
seed poloidal field. For the hydrodynamic SSP, the three- 
dimensional instabilities are inflectional instabilities of the 
spanwise-modulated finite-amplitude streaks. In the Ke- 
plerian MHD problem, a natural instability candidate is a 
three-dimensional MRI of the (9(1) toroidal magnetic field. 

We use a three-dimensional nonlinear continuation code 
based on Newton iteration to test this scenario quantita- 
tively. The code is similar to those used to study hydrody- 
namic SSPs lflol,IT2ll and has been tested extensively with 
standard nonlinear problems HjIi . Collocation on a Gauss- 
Lobatto grid is used in y and a Fourier representation (with 
dealiasing) is used in x and z- The Newton solver relies on 
the linear algebra library LAPACK to solve real systems 
with (9(20 000) unknowns. An iterative generalized eigen- 
value problem solver from the ARPACK library is used 



to address 2D linear stability problems. The convergence 
of our solutions (~ 10"'^ for energy) has been thoroughly 
checked by looking at Chebyshev and Fourier spectra and 
by comparing results obtained at diff'erent resolutions. 

We adopt a forcing strategy ifioll . which first consists in 
artificially forcing nonlinear steady 2D-3C solutions and in 
computing their three-dimensional instability modes. One 
must then check that the nonlinear feedback of these modes 
can take over the forcing to obtain three-dimensional non- 
linear unforced solutions. For the MHD problem, an ap- 
propriate forcing is a toroidal electromotive force (EMF) 



EMF,(3;,z) = 
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i|^^jcos/?z. 
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applied to Eq. Q with A = (9(1), which compensates for 
the resistive decay of an (9(1 /Rm) poloidal magnetic field 
that in turn generates an (9(1) toroidal field. Fig. [T] depicts 
a 2D-3C magnetic field solution of the forced nonlinear 
MHD equations obtained with our Newton solver 




FIG. 1: 2D-3C magnetic configuration induced by the artificial 
toroidal EMF ((S}, /? = 1 and perfectly conducting boundary con- 
ditions, bx is represented on a color scale (black to white from - 
max to max) and (by, b-) by arrows. Rm = 750, Re = 10, A = 1.5, 
max(Z7,.) = 1.006 and max(fo,) = 0.00195. (Ny, N^) = (32, 32). 

We then consider the stability of forced 2D-3C MHD 
flows with a dominant (9(1) toroidal magnetic field to in- 
finitesimal perturbations with exp(/Q'x) dependence, and 
compute the associated MRI eigenmodes with our linear 
eigenvalue solver (the Alfven continuum is removed here 
owing to viscous and resistive effects). For the forcing 
term (O steady solutions have three symmetries which 
are used to reduce computational costs: first, they either 
have reflect (R) z — > -z or shift-and-reflect (SR) (x, z) — > 
(x + Lx/2, -z) symmetry which both turn [v^, Vy, v.) into 
(v^, Vy, -Vj,) and (b^, by, b.) into {-b^, -by, bj. They also 
have either z-shift (S) z -» z H- LJ2 or double-shift (DS) 
symmetry (x, z) — » (x -H Lv/2,z + L-/2) which both turn 
(v, b) into (v, -b). Finally, owing to the invariance of 
PCF under z-rotations by n, all solutions have shift-and- 
rotate (x, y, z) — » iLx/2 - x, -y, z + L^I2) symmetry, which 
turns fields (/i v, hy, h~) into (-/i v, -hy, h,) (all transforma- 
tions require appropriate x and z phase choices). The a- 
dependence of the growth rates of several modes for the 
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2D-3C configuration of Fig. [T]is plotted in Fig. |2] A (y, z) 
cut through the marginal (R,S) mode is shown in Fig. [3] 



£ 



S 0, 



C 0.01 H 



0.00 




0.0 



0.1 0.2 0.3 0.4 

Streamwise wavenumber a 



FIG. 2: Growth rates of the most unstable three-dimensional MRI 
eigenmodes for the 2D-3C configuration with 0(1) toroidal field 
shown in Fig.[T] as a function of the streamwise wavenumber a. 
The legend indicates the symmetries of each mode. 
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FIG. 3: Magnetic field in the x = plane of the near marginal 
(R,S) MRI eigenmode with a = 0.375, depicted by a black dot in 
Fig. [2] Same representation and parameters as in Fig.[T] 

In order to find a good initial guess to discover steady 
three-dimensional solutions with the Newton solver, one 
must select a marginally stable mode and check whether 
or not the jc-averaged EMF due to its nonlinear interac- 
tions can take on the role of the artificial EMF (O to sus- 
tain the poloidal field. This depends strongly on the se- 
lected MRI mode and on its streamwise wavenumber a. 



much like in the hydrodynamic problem 111211 . If the feed- 
back is bad, the amplitude A of the artificial EMF (and 
therefore the toroidal component of the 2D-3C field) is 
adjusted and a 2D-3C flow marginally unstable to a MRI 
mode with slightly different a is recalculated. This way, 
we identify a (R,S) MRI eigenmode (Fig. [3]) whose non- 
linear interactions create a toroidal EMF that has an inter- 
esting positive correlation with the artificial EMF (Fig.lH. 
Note that considering the feedback in the 2D-3C momen- 
tum equation is not important to find an interesting MRI 
mode since the goal is to replace an artificial forcing term 
imposed in the 2D-3C induction equation only. We final at- 
tempt a continuation with respect to the forcing amplitude 
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FIG. 4: Full line: normalized shearwise profile of the k- = fH 
component of the toroidal EMF created by the self-interactions 
of the three-dimensional a = 0.375 MRI mode of Fig.|3] Dashed 
line: normalized shearwise profile of the artificial EMF l[3j. 



to find unforced (A = 0) solutions. The code is initialized 
with the 2D-3C base flow of Fig. [T]plus a small amount 
of the three-dimensional (R,S) MRI mode of Fig.|3] and A 
is set to a slightly smaller value than that used to force 
the 2D-3C flow. The solver converges to a fully three- 
dimensional forced solution, demonstrating that the bifur- 
cation to three-dimensional solutions is subcritical with 
respect to A. Further confirmation of the crucial role of 
the MRI mode feedback is obtained by performing simi- 
lar calculations initialized with marginal MRI modes with 
negative feedback (such as the (SR,DS) mode of Fig. |2] 
at Q- = 0.355), which reveal supercritical behaviour. The 
subcritical (R,S) branch can be followed down to A = 
(Fig. Is), at which point (black dot) nonlinear interactions 
due to the x-dependent part of the solution fully take over 
the forcing term. The A = point is therefore a fully 
three-dimensional subcritical nonlinear steady solution of 
the original MHD equations with no forcing (Fig.|6]l. The 
x-dependence of the solution, unlike for A = 1.5, can no 
longer be described by a single Fourier mode. 
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FIG. 5: y-integrated amplitude of the = 0.375, k- = com- 
ponent of v,, for three-dimensional steady forced solutions, as a 
function of the forcing amplitude A. (N_„Ny,N^) = (12,32,32). 
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FIG. 6: Cuts through the unforced (A = 0) nonlinear steady so- 
lution of Fig. [5] (a = 0.375, black dot). Top: h, a.t z = L,/2. 
Bottom: fe, (color scale) and {b^, b-) (arrows) at x = L^IA. 
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FIG. 7: Continuation with respect to Rm for Re = 10, (a,P) = 
(0.375, 1). Inset: continuation with respect to Re for Rm = 1500 
(full Une), Rm = 3000 (dashed line). (A'.,, A',,, M) = (8, 24, 32). 



As shown in Fig. |2l this nonlinear steady solution can 
be continued to large Rm but exists only in a narrow range 
of Re. The reason for this could be that only at low Re do 
marginally unstable MRI modes of the toroidal magnetic 
field originate from a steady bifurcation. At larger Re, 
marginal MRI modes systematically arise from a Hopf bi- 
furcation instead: pairs of steady eigenvalues correspond- 
ing to modes with the same symmetry collide and turn into 
complex-conjugate pairs for increasing Re. To discover 
similar dynamo solutions at larger Re, it may therefore 
prove necessary to consider time-dependent MRI modes 
instead of steady ones. This is unfortunately numerically 



far more challenging, since the symmetries used to reach 
decent three-dimensional resolutions are broken when the 
MRI modes turn into complex-conjugate pairs. 

We have presented an instance of self-sustaining, non- 
linear dynamo solution in Keplerian PCF, whose critical 
Rm ^ 670 is comparable to that found in zero-net-flux 
simulations (another Rm definition is used in iQl). Prelim- 
inary direct spectral numerical simulations seem to con- 
firm independently the existence of this solution, whose 
structure is dominated by a coherent zero-net-flux toroidal 
magnetic field. We also discovered a (R,DS) branch, which 
makes it probable that many solutions with different sym- 
metries exist. Such coherent structures are strictly speak- 
ing not turbulence but, like the hydrodynamics SSP, they 
probably act as organizing centers of the dynamics in phase 
space ilOr and could play an important role in triggering 
and sustaining MHD turbulence in magnetized Keplerian 
disks, where all the basic ingredients of the dynamo are 
present. Hopefully, this idealized study will be helpful to 
uncover the details of the dynamo in more realistic set-ups. 
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